Jet suppression of pions and single electrons at Au+Au collisions at RHIC 
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Jet suppression is considered to be a powerful tool to study the properties of a QCD medium 
created in ultra-relativistic heavy ion collisions. However, theoretical predictions obtained by using 
jet energy loss in static QCD medium show disagreement with experimental data, which is known 
as the heavy flavor puzzle at RHIC. We calculate the suppression patterns of pions and single 
electrons for Au+Au collisions at RHIC by including the energy loss in a finite size dynamical QCD 
medium, with finite magnetic mass effects taken into account. In contrast to the static case, we here 
report a good agreement with the experimental results, where this agreement is robust with respect 
to magnetic mass values. Therefore, the inclusion of dynamical QCD medium effects provides a 
reasonable explanation of the heavy flavor puzzle at RHIC. 

PACS numbers: 12.38.Mh; 24.85.+p; 25.75.-q 
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INTRODUCTION 

Jet suppression [l[ measurements at RHIC and LHC, 
and their comparison with theoretical predictions, pro- 
vide a powerful tool for mapping the properties of a 
QCD medium created in ultra-relativistic heavy ion col- 
lisions 0-3- However, jet suppression predictions, done 
under assumption of static QCD medium showed a 

disagreement with the available data from RHIC experi- 
ments 

mm 

. This disagreement has been named "heavy 
flavor puzzle at RHIC" [TH, EH, and raised important 
questions about the ability of the available theories to 
model the matter created at ultra-relativistic heavy ion 
collisions at RHIC. 

Since the suppression results from the energy loss of 
hig h energy partons moving through the plasma 
Il7| . accurate computations of jet energy loss mechanisms 
are essential for the reliable predictions of jet suppres- 
sion. In [l8|,[lj|, we developed a theoretical formalism for 
the calculation of radiative energy loss in realistic finite 
size dynamical QCD medium (sec also a viewpoint [Hj]); 
which abolished a static approximation used in previous 
models j20N26ll . Furthermore, in 27j], we extended the 
study from [18[ to include a possibility for existence of 
finite magnetic mass; this generalization was motivated 
by various non-perturbative approaches [28 - 31 [ , which 
report non-zero magnetic mass. These studies, together 
with the previously developed collisional energy loss for- 
malism in finite size dynamical QCD medium [32(, enable 
us to provide the most reliable computations of the en- 
ergy loss in QGP so far. 

In this paper, we integrate the developed energy loss 
formalism into a computational framework that can gen- 
erate reliable predictions for RHIC and LHC experimen- 
tal data. The numerical procedure includes: i) both colli- 
sional and radiative energy loss from the n ewly developed 
(dynamical QCD medium) formalism [H [ll E3, H3], ii) 
multi-gluon fluctuations, i.e. the fact that energy loss 
is a distribution [33j . and iii) path length fluctuations, 



i.e. the fact that particles travel different paths in the 
medium [34j |. We use this framework to generate sup- 
pression predictions for pions and single electrons at most 
central 200 GeV Au+Au collisions at RHIC. The gener- 
ated predictions are directly compared with RHIC exper- 
imental data in order to test our understanding 
of QGP created at these collisions. 



COMPUTATIONAL FRAMEWORK 



The quenched spectra of partons, hadrons, and lep 



tons are calculated as in 134 
convolution 



35j from the generic pQCD 



E f d 3 a(e) Eid s a(Q) 



dp 3 . 



dpi 



P(Ei -> E f ) 
Hq) ® f(H Q -> e), 



(1) 



where Q denotes quarks and gluons. For charm and bot- 
tom, the initial quark spectrum, Eid 3 cr(Q)/dp 3 ^ is com- 
muted at ncxt-to-leading order using the code from [36|, 
37| : for gluons and light quarks, the initial distributions 
are computed at leading order as in [lH • P{E% — > Ef) 
is the energy loss probability, D(Q — > Hq) is the frag- 
mentation function of quark or gluon Q to hadron Hq. 
The last step (/(Hq — > e)) is only applicable for heavy 
quarks, and it represents the decay function of hadron 
Hq into the observed single electron. We use the same 
mass and factorization scales as in [39| and employ the 
CTEQ5M parton densities with no intrinsic kx- As 
in [39l | we neglect shadowing of the nuclear parton distri- 
bution. 

We assume that the final quenched energy Ef is large 
enough that the Eikonal approximation can be employed. 
We also assume that in Au+Au collisions, the jet to 
hadron fragmentation functions are the same as in e + e~ 
collisions. This assumption is expected to be valid in 
the deconfined medium case, where hadronization of 
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Q — > Hq cannot occur until the quark emerges from the 
QGP. 

As in [34|], the energy loss probability P(Ei — > Ef) 
is generalized to include both radiative and collisional 
energy loss and their fluctuations. However, a major dif- 
ference between [3~i ] and the present study is that we 
here take into account both the radiative \lm and colli- 



sional 32j energy losses in a realistic finite size dynamical 
QCD medium. 

To take into account geometric path length fluctua- 
tions in the energy loss probability, we use [34|: 

P( Ei ->E f = Ei - A rad - A coll ) = 

J dh P(L) P rad {A rad] L) ® P coll (A coll ; L) . (2) 

Here P(L) is the distribution of path lengths traversed by 
hard scatterers in 0-5% most central collisions, in which 
the lengths are weighted by the probability of production 
and averaged over azimuth. Note that currently two def- 
initions for P(L) (see [Ilj]), one from [34| and the other 
from [42j] , are commonly used. Since these distributions 
are significantly different (see [4l|). we will use both of 
them in the analysis. Also, since P(L) is a purely geo- 
metric quantity, it is the same for all jet varieties. 

P rad (A rad ;L) and P coU {A coU ; L) in Eq. © are, re- 
spectively, the radiative and collisional energy loss prob- 
abilities. The procedure for including fluctuations of 
the radiative energy loss probability (P r ad(A rad ] L)) due 
to gluon number fluctuations is discussed in detail in 
Ref. [Ill 43 1 . Note that the procedure is here gener- 
alized to include the radiative energy loss in finite size 
dynamical QCD medium 18, [Hj], as well as a possibil- 
ity for existence of finite magnetic mass 27[; in par- 
ticular, we extract the gluon radiation spectrum from 
Eq. (fO) in [27|. For collisional energy loss probability 
{Pcoii(A co u; L)), the full fluctuation spectrum is approxi- 
mated by a Gaussian centered at the average energy loss 

coll ( 



with variance a 2 coll = 2T(AE coU (pj_, 
AE coll (p±,L) is given by Eq. (14) in H 



Here 



T is the tem- 
perature of the medium, p±^ is the initial momentum of 
the jet, and L is the length of the medium traversed by 
the jet. 

We note that, in the suppression calculations, we sepa- 
rately treat radiative from collisional energy loss; Conse- 
quently, we first calculate the modification of the quark 
and gluon spectrum due to radiative energy loss, and 
then due to collisional energy loss in QCD medium. This 
is a reasonable approximation when the radiative and 
collisional energy losses can be considered small (which 
is in the essence of the soft-gluon, soft-rescattering ap- 
proximation used in all energy loss calculations sofar [18l — 
261 |45|), and when collisional and radiative energy loss 
processes are decoupled from each other (which is the 
case in the HTL a ppro ach (4(| used in our energy loss 
calculations [H 0, 152] ) . Also , we assume that strong 
coupling constant as is fixed at 0.3. 



Finally, to obtain 7r° suppression from quark and gluon 
suppression, we use the following estimate [13, [48| 

R AA (ir°, P± ) a f g R AA (g,p±) + (1 - f g )R AA (l, Pl _), (3) 



where f„ 



-pj_/10.5GeV 



is the fraction of pions with a 



given momentum p± that arise from gluon jet fragmenta- 
tion, R AA (g,p±) is the gluon suppression and R AA (l,p±) 
is the light quark suppression. 

To calculate D and B meson suppression, we will use 
both delta and Peterson fragmentation functions. 
Furthermore, to obtain single electron suppression, we 
will use the following estimate [||: 



RAA(e ± ,P±) ~ 
R AA (D,2p ± )^EA 



dpi 



R AA (B,2 P± ) 



da B (2p 1 



dcr D (2p ± ) . d<y B (2p±_) 



(4) 



where da d p^ i ^dp^ ) * s ^ (^) meson initial momen- 
tum distribution, and R AA (D,p±) (R AA (B,p±)) is the 
D (B) meson suppression. 

NUMERICAL RESULTS 

In this section, we concentrate at 200 GeV Au+Au 
collisions at RHIC, and present our suppression pre- 
dictions for light and heavy flavor observables. For- 
tius, we consider a quark-gluon plasma of temperature 
T = 225MeV, with iV/=2.5 effective light quark flavors 
and strong interaction strength as = 0.3, as representa- 
tive of average conditions encountered in Au+Au col- 
lisions at RHIC. For the light quarks we assume that 
their mass is dominated by the thermal mass 
where [i = gT ^Jl+N j /Q « 0.5 GeV is the Debye screen- 
ing mass. The gluon mass is taken to be m g = 
I1/V2. For the charm (bottom) mass we use M = 1.2 GeV 
(M = 4.75 GeV). 

Figure [T] shows momentum dependence of quarks and 
gluon suppressions at RHIC, obtained by using path 
length distributions from j34|. We observe a clear hi- 
erarchy between the quarks and gluon suppressions: i) 
bottom quark is significantly less suppressed than charm 
quark; ii) charm and light quarks have similar suppres- 
sions for initial jet energies larger than 5 GeV; Hi) glu- 
ons are significantly more suppressed than all types of 
quarks. This already observed/established hierarchy (see 
e.g. [43|) therefore remains valid for the case of dynam- 
ical QCD medium as well, despite the fact that both 
inclusion ofpath length fluctuations (Hf and dynami- 
cal effects [19( into the suppression calculations tend to 
reduce the difference between different types of quarks 
and gluons. We also observe that inclusion of magnetic 
mass can decrease the jet suppression (for all types of 
quarks and gluons) from 25-50%, compared to the case 
of zero magnetic mass. While in Fig. [1] we show only the 
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FIG. 1: Quarks and gluon suppressions are presented as a function of initial jet energy for 200GeV Au+Au collisions at RHIC. 
The panels are obtained by using path length distributions from [33]. On each panel full curves correspond to the case when 
magnetic mass is equal to zero. Gray bands correspond to the case when magnetic mass is non-zero (i.e. 0.4 < hm/he < 0.6 [U 
EH!), where the lower boundary corresponds to [im/he = 0.4 and the upper boundary corresponds to /im/^b = 0.6. 
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Left panels show the compari- 
pion suppression predictions with 
[Tol ] experimen- 



FIG. 2: 

son of 

7T° PHENIX g| and STAR 
tal data from 200 GeV Au+Au collisions 
at RHIC. Right panels show the compari- 
son of single electron suppression predictions 
with non-photonic single electron data from 
PHENIX [| and STAR [ll[ at 200 GeV 
Au+Au collisions. For the two upper pan- 
els, and the two lower panels, the suppres- 
sion predictions are obtained, respectively, 
by using path length distributions from [34| 
(marked with "W") and [12] (marked with 
"D"). On each panel, the gray region cor- 
responds to the case when (xm > (i.e. 
< fiAi/fiE < 0.6), where the lower bound- 
ary corresponds to hm / (Xe = and the up- 
per boundary corresponds to hm/he = 0.6. 



results by using path length distributions from 34j . we 
observe somewhat lower suppression results when distri- 
bution from 42j is used (data not shown); this is expected 
ha ving in mind higher probability for lower path-lengths 



421 compared to [34| (see |4l|) 



To calculate D and B meson suppression, we used both 
delta and Peterson [4o| ] fragmentation functions, and ob- 
served that the choice of fragmentation function only 
marginally changes the value of suppression Raa (data 
not shown). The reason behind the small difference is 
that fragmentation functions do not significantly mod- 
ify the distribution slopes [lU, due to which the meson 
suppression becomes insensitive to the choice of fragmen- 
tation function. For simplicity, we will therefore further 
use delta function fragmentation for the calculation of 
single electron suppression. 

Figure [2] shows momentum dependence of pion and sin- 
gle electron suppressions at RHIC, obtained by using two 
path length distributions (from [34| and from 42J). The 
predictions are compared with the relevant PHENIX 0, 
[U and RHIC 1(J U| experimental data at Au+Au col- 
lisions at RHIC. For both path length distributions we 



observe a reasonable agreement with the experimental 
data, which is moreover robust with respect to non-zero 
magnetic mass introduction. 

Furthermore, from Fig. [2] we see that at jet energies 
above 15 GeV, the predicted pion and single electron 
suppressions become very similar. This prediction is rea- 
sonable, since i) at high jet energies suppression patterns 
for all types of quarks become similar, and ii) above 
10 GeV pion suppression is strongly dominated by light 
quark suppression (i.e. gluon contribution to pion sup- 
pression becomes negligible) [13, |H|. This behavior is 
qualitatively different from the one below 10 GeV, where 
obtained suppressions are notably different (by < 2). It 
will be interesting to compare this predicted pattern with 
the upcoming high luminosity RHIC data. 



CONCLUSIONS 

In this paper we calculated the suppression pattern 
of pions, D and B mesons and single electrons in cen- 
tral 200 GeV Au+ Au collisions at RHIC energies. The 
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calculation is based on the radiative and collisional en- 
ergy loss in a finite size dynamical QCD medium, which 
is a key ingredient for obtaining reliable predictions for 
jet quenching in ultra-rclativistic heavy ion collisions. 
This energy loss formalism was here integrated into a 
computational framework that includes multi-gluon and 
path length fluctuations. We obtained a reasonably good 
agreement between the generated suppression patterns 
and experimental data at Au+Au collisions at RHIC, and 
this agreement is robust with respect to introduction of 
finite magnetic mass. The agreement strongly suggests 
that the main deficiency responsible for the "heavy fla- 
vor puzzle" at RHIC was the static approximation, i.e. 
the fact that dynamical nature of plasma constituents 
was not taken into account. Predictions of the dynami- 
cal energy loss formalism remain to be tested against the 
upcoming high luminosity RHIC and LHC data. 

Acknowledgments: This work is supported by Marie 
Curie International Reintegration Grant within the 7 th 
European Community Framework Programme (PIRG08- 
GA-2010-276913) and by the Ministry of Science and 
Technological Development of the Republic of Serbia, un- 
der projects No. ON171004 and ON173052. 



[1 

[2 
[3 
[4 

[5 

[6 

[7: 

[8 
[9 
[10 

[11 

[12 
[13 
[14 

[15 
[16 
[17 



J.D. Bjorken: FERMILAB-PUB-82-059-THY (1982). 

N. Brambilla et aL, Preprint [ hep-ph/0412158] (2004). 

M. Gyulassy, Lect. Notes Phys. 583, 37 (2002). 

D. d'Enterria, B. Betz, Lect. Notes Phys. 785, 285 

(2010). 

M. Djordjevic, M. Gyulassy, S. Wicks and R. Vogt, Phys. 
Lett. B 632, 81 (2006). 

R. Rapp, V. Greco and H. van Hees, Nucl. Phys. A 774, 
685 (2006). 

N. Armesto, A. Dainese, C. A. Salgado and U. A. Wiede- 
mann, Phys. Rev. D 71, 054027 (2005). 
A. Adare et al. [PHENIX collaboration], Phys. Rev. Lett. 
101, 232301 (2008). 

A. Adare et al. [PHENIX collaboration], Phys. Rev. Lett. 
98, 172301 (2007). 

B. I. Abelev et al. [STAR collaboration], Phys. Rev. C 80 
44905 (2009). 

B.I. Abelev et al. [STAR collaboration], Phys. Rev. Lett. 
98, 192301 (2007). 
M. Gyulassy, Physics 2, 107 (2009). 
M. Djordjevic, J. Phys. G 32, S333 (2006). 
M. Guylassy, I. Vitev, X. N. Wang and B. W. Zhang, in 
Quark Gluon Plasma 3, edited by R. C. Hwa and X. N. 
Wang, p. 123 (World Scientific, Singapore, 2003) 
R. Baier, Yu. L. Dokshitzer, A. J. Mueller and D. Schiff, 
Phys. Rev. C 58, 1706 (1998). 
R. Baier, D. Schiff, B. G. Zakharov, Ann. Rev 
Part. Sci. 50, 37 (2000). 

A. Kovner and U. A. Wiedemann, in Quark 
Plasma 3, edited by R.C. Hwa and X.N. Wang, 



Nucl. 

Gluon 
p. 192 



[18 
[19 

[20 

[21 

[22 

[23 
[24 

[25 

[20 
[27 
[28 

[29 

[30 

[31 

[32 
[33 

[34 

[35 

[36 

[37 

[38 

[39 
[40 

[41 

[42 
[43 

[44 
[45 



[46 
[47 
[48 
[49 



(World Scientific, Singapore, 2003). 

M. Djordjevic, Phys. Rev. C 80, 064909 (2009). 

M. Djordjevic and U. Heinz, Phys. Rev. Lett. 101, 

022302 (2008). 

M. Djordjevic and M. Gyulassy, Phys. Lett. B 560, 37 
(2003); and Nucl. Phys. A 733, 265 (2004). 
M. Gyulassy, P. Levai and I. Vitev, Nucl. Phys. B 594, 
371 (2001). 

M. Gyulassy and X. N. Wang, Nucl. Phys. B 420, 583 
(1994); X. N. Wang, M. Gyulassy and M. Plumer, Phys. 
Rev. D 51 (1995) 3436. 

U.A. Wiedemann, Nucl. Phys. B 588, 303 (2000); and 
Nucl. Phys. B 582, 409 (2000). 

E. Wang and X. N. Wang, Phys. Rev. Lett. 87 142301, 

(2001) . X. N. Wang and X. F. Guo, Nucl. Phys. A 696, 
788 (2001); X. F. Guo and X. N. Wang, Phys. Rev. Lett. 
85 3591, (2000). 

N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. 
Rev. D 69, 114003 (2004). 

M. Djordjevic, Phys. Rev. C 73 , 044912 (2006). 
M. Djordjevic, larXiv:1105.4359l [nucl-th] 
Yu. Maezawa et al. [WHOT-QCD Collaboration], Phys. 
Rev. D 81 091501 (2010); Yu. Maezawa et al. [WHOT- 
QCD Collaboration], PoS Lattice 194 (2008). 
A. Nakamura, T. Saito and S. Sakai, Phys. Rev. D 69, 
014506 (2004). 

A. Hart, M. Laine and O. Philipsen, Nucl. Phys. B 586, 
443 (2000). 

D. Bak, A. Karch and L. G. Yaffe, JHEP 0708, 049 
(2007). 

M. Djordjevic, Phys. Rev. C 74, 064907 (2006). 

M. Gyulassy, P. Levai and I. Vitev, Phys. Lett. B 538, 

282 (2002). 

S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, 
Nucl. Phys. A 784, 426 (2007). 

M. Djordjevic, M. Gyulassy, R. Vogt and S. Wicks, Phys. 
Lett. B 632 81 (2006). 

M. Cacciari, P. Nason and R. Vogt, Phys. Rev. Lett. 95 
122001 (2005). 

M. L. Mangano, P. Nason and GRidolfi, Nucl. Phys. B 
373, 295 (1992). 

I. Vitev and M. Gyulassy, Phys. Rev. Lett. 89, 252301 

(2002) . 

R. Vogt, Int. J. Mod. Phys. E 12, 211 (2003). 

H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 
12, 375 (2000). 

S. Wicks, PhD thesis, AAT-3333486, PROQUEST- 

1614268891, 2008. 338pp 

A. Dainese, Eur. Phys. J. C 33 495 (2004). 

M. Djordjevic, M. Gyulassy and S. Wicks, Phys. Rev. 

Lett. 94, 112301 (2005). 

G. D. Moore, D. Teaney, Phys. Rev. C 71, 064904 (2005). 
P. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0111, 
057 (2001); JHEP 0206, 030 (2002); JHEP 0301, 030 

(2003) . 

M. Djordjevic, Nucl. Phys. A 783 197 (2007). 

I. Vitev, Phys. Lett. B 606, 303 (2005). 

A. Adil and M. Gyulassy, Phys. Lett. B 602 (2004) 52. 
C. Peterson, D. Schlatter, I. Schmitt and P. M. Zerwas 
Phys. Rev. D 27, 105 (1983). 



